3D Lattice-Boltzmann Model for Magnetic Reconnection 
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In this paper we develop a 3D Lattice-Boltzmann model that recovers in the contirmous limit the 
two-fluids theory for plasmas, and consecuently includes the generalizated Ohm's law. The model 
reproduces the magnetic reconnection process just by given the right initial equlibrium conditions 
in the magnetotail, without any assumption on the resistivity in the diffusive region. In this model, 
the plasma is handled like two fluids with an interaction term, each one with distribution functions 
associated to a cubic lattice with 19 velocities (D3Q19). The electromagnetic fields are considered 
like a third fluid with an external force on a cubic lattice with 13 velocities (D3Q13). The model 
can simulate either viscous fluids in the incompressible limit or non- viscous compressible fluids, and 
sucessfuUy reproduces both the Hartmann flow and the magnetic reconnection in the magnetotail. 
The reconnection rate obtained with this model is -R=0.109, which is in excellent agreement with 
the observations. 

PACS numbers: 94.30.cp, 52.30.Ex, 52.65.-y 
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I. INTRODUCTION 

The magnetic reconnection is one of the most interest- 
ing phenomenon of plasma physics. This process quickly 
transforms the magnetic energy into termic and kinetic 
energies of the plasma. It is mostly observed inside 
of astrophysical plasmas, such as solar flares (where it 
contributes to the plasma heating), and in the terres- 
trial magnetosphere, where it support the income flux of 
plasma and electromagnetic energy. 

The magnetic reconnection requires the existence of a 
diffusive region, where dissipative electric fields change 
the magnetic field topology. The first models were inde- 
pendently formulated by Sweet [J, in 1958, and Parker 
[2I, in 1957. They suggested that the magnetic recon- 
nection is a steady-state resistive process that occurs in 
the vicinity of a neutral line. This model reduces the 
phenomenon to a boundary condition problem and can 
explain the magnetic field reconnection. However, it has 
some problems when compared with experimental obser- 
vations (i.e. a very slow reconnection rate), and it leaves 
unexplained the origin of the high-resistive region. In 
1964, Petschek Q proposed the first model for fast re- 
connection rates. He included a much smaller diffusion 
region than the Sweet-Parker model, but he suggested 
that the rest of the boundary layer region should consist 
of slow shock waves that accelerate the plasma up to the 
Alfven velocity. Nevertheless, the origin of the diffusive 
region remains unexplained. 
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At present, the nature of this phenomenon has been 
studying by using kinetic theory and considering coUi- 
sionless plasmas, since this is a common property of as- 
trophysical plasmas. One of the developments of the ki- 
netic theory is the generalized Ohm's law, where some 
extra terms explain the existence of a dissipative electric 
field. The introduction of these extra terms in resistive 
magnetohydrodynamics is called MHD-Hall jj] . A useful 
approximation of the kinetic theory consists of modelling 
the plasma like two fluids (one electronic and one ionic), 
which have independent momentum, mass conservation 
and state equations, plus an interaction term in the mo- 
mentum equation |j]. This treatment, in the one- fluid 
limit, introduces in a natural way the extra terms of the 
generalized Ohm's law. However, the equations involved 
by this treatment are complex and it is difficult to find 
an analytic solution for any problem. 

For this reason, most plasma processes are studied by 
numerical methods. One of the numerical methods for 
simulating fiuids is Lattice Boltzmann (LB) 5], which 
was developed from lattice-gas automata. Lattice Boltz- 
mann simulations are performed on regular grids of many 
cells and a small number of velocity vectors per cell, each 
one associated to a density distribution function, which 
evolve and spread together to the neighbohr cells accord- 
ing to the coUisional Boltzmann equation. The first LB 
model for studying plasmas reproduces the resistive mag- 
netohydrodynamic equations and was developed by Chen 
[a, Q as an extension of the Lattice-Gas model devel- 
oped by Chen and Matthaeuspj and Chen, Matthaeus 
and Klein l!||. This LB model uses 37 velocity vectors 
per cell on a square lattice and is developed for two di- 
mensions. Thereafter, Martinez, Chen and Matthaeuss 
|10l | decreased the number of velocity vectors from 37 



to 13, which made easier a future 3D extension. One 
of the first LB models for magnetohydrodynamics in 3D 
was developed by Bryan R. Osborn in his master thesis 
[ll|. He used 19 vectors on a cubic lattice for the fluid, 
plus 7 vectors for the magnetic field, which makes a to- 
tal number of 26 vectors per cell. By following a different 
path, Fogaccia, Benzi and Romanelli 12] introduced a 3D 
LB model for simulating turbulent plasmas in the elec- 
trostatic limit. All these models reproduce the resistive 
magnetohydrodynamc equations for a single fluid. 

In this paper, we introduce a 3D Lattice-Boltzmann 
model that recovers the plasma equations in the two- 
fluids theory. In this way, the model is able to reproduce 
magnetic reconnection, without the a priori introduction 
of a resistive region. Moreover, it is able to reproduce the 
fluid state-equation with a general polytropic coefficient. 
The model uses 39 vectors per cell and 63 probability 
density functions (19 for each fiuid, 25 for the electrical 
and magnetic fields). In section^ we describe the model, 
with the evolution rules and the equilibrium expressions 
involved for the 63 density functions, plus the way to 
compute the electric, magnetic and velocity fields. The 
Chapman-Enskog expansion showing how these rules re- 
cover the two-fiuids magnetohydrodynamic equations is 
developed in Appendix^ In order to validate the model, 
we simulate the 2D Hartmann's fiow in section ITTll and, 
finally, the magnetic reconnection for a magnetotail equi- 
librium configuration in section llVl The main results and 
conclusions are summarized in section Ivl 



II. 3D LATTICE-BOLTZMANN MODEL FOR A 
TWO-FLUIDS PLASMA 

In a simple Lattice-Boltzmann model [5|, the D- 
dimensional space is divided into a regular grid of cells. 
Each cell has Q vectors Vi that links itself with its neigh- 
bors, and each vector is associated to a distribution func- 
tion fi. The distribution function evolves at time steps 
6t according to the Boltzmann equation. 




FIG. 1: Cubic Lattice D3Q19 for modelling the electronic and 
ionic fluids. The arrows represent the velocity vectors vf and 
p indicates the plane of location. 




Mx + v,6t, t + 6t)^ f,{x, t) = a,(f, t) 



(1) 



FIG. 2: Cubic Lattice D3Q13 for modelling the electric field. 
The arrows represent the electric vectors e*" . 



electric field (figurelSJ and 7 vectors for the magnetic field 
(figure |3J- The velocity vectors are denoted by wf, where 
i = 1,2,3,4,5,6 indicates the direction and p — 0,1,2 
indicates the plane of location. Their components are 

v^^ = cV2{cos{{2i - l)7r/4), sin((2i - l)7r/4), 0) , (3a) 



where fli {x, t) is a collision term, which is usually taken 
as a time relaxation to some equilibrium density, f'jl'^. 
This is known as the the Bhatnagar-Gross-Krook (BGK) 
operator p!a |. 



n,{x,t)^^-{Mx,t)^fP{x,t)) 

T 
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where r is the relaxation time and f°'^{x, t) is the equilib- 
rium function. The equilibrium function is chosen in such 
a way, that (in the continuum limit) the model simulates 
the actual physics of the system. 

For our 3D model, we use a cubic regular grid, with 
lattice constant 5x=^/2c5t and c is the light speed 
(c~3 X lO^m/s). There are 19 velocity vectors for the 
electronic and ionic fluids (figure ^1, 13 vectors for the 



v] = c\/2(cos((2i - l)7r/4), 0, sin((2i - l)7r/4)) , (3b) 

vf = c%/2(0, cos((2i - l)7r/4), sin((2i - l)7r/4)) , (3c) 
for i < 5, and 

= cV2((-l)\0,0) , (4a) 
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/i=cV2(0,(-l)\0) 



i/2=cV2(0,0,(-l)^) 



(4b) 



(4c) 




as follows: 



FIG. 3: Cubic Lattice D3Q7 for simulating the magnetic field, 
the arrows indicate the magnetic vectors 6f . 
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FIG. 4: Index relationship between the velocity vectors and 
the electric and magnetic vectors. 



for i > 5. This makes 18 vectors. The missing one is the 
rest vector v^, with componets (0,0,0). 

The set of 13 electric field vectors, ef-, and 7 magnetic 

field vectors, h\, are related with the velocity vectors as 
follows: 
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where the index i takes the values i=l, 2, 3, 4. 

The distribution functions that describe the fluids, de- 
noted by /f and /g , propagate with each velocity vec- 
tor v^ and with the rest vector vq, respectively, and uses 
these vectors to compute the velocity fields for each fluid. 
Here, the index s distinguishs between electronic (s=0) 
and ionic {s—\) fluids. Similarly, the distribution func- 
tions associated for the electromagnetic field are denoted 
by /j^- and /q . They also propagate in the direction 
of the velocity vectors v^ and zTq, but they use the elec- 
tric and magnetic field vectors to compute those fields. 
Summarizing, The macroscopic variables are computed 
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where ps and Vg are the density and velocity of each fluid, 
and rris and q^ are its particle mass and charge (here, s=0 
represents electrons and s=l represents ions, as before). 
In addition, E and B are the electric and magnetic fields, 
J is the total current density and pc is the total charge 
density. 

For their evolution, we follow the proponsal of J.M. 
Buick and C.A. Created for the lattice Boltzmann equa- 
tions [3, 
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i''\x,t + 6t)-i''\x,t) = n[''\x,t) , (10) 
where K — 0,1, 2. The force vectors F^*^ in Eq.® are 
F'^'^^^PsiE + Vs^B) 



is) 



lypsiVs - V(^s+l)mod2) + ^0 ' 



(11) 



is) 



where v is the collision frequency of the plasma, Fq 
is any external force (for instance, a gravitational force) 
and the equilibrium density current vector J' in Eq. 
is defined by 






Ps 



(12) 



The collision terms il^- ' and fij, ' are given by 



-)P(s) 



fPis)/ 



rP{s)cq, 



^r' = --ifn^^t) - fr^i^^t)) , (13a) 



incompressible and viscous fluids. The first ones are gov- 
erned by the continuity equation 



v-(p.K') + ^-o , 



(19) 



-,P(2) 



^r = --(/r (^^'*) - nr i^^t)) > (ist) 



^f^ = --(/r(-^,t)-#>'^(f,0) , (13c) 

where tk is the relaxation time, kk= ?^7 and As = ^ 



2rA' 



The equilibrium functions for the fluids, f[ and 



,p(s)oq 



2t^, ■ 



/(i^^^'^ are 



j,p(s)oq 



/r'^^'^(f,i)=c.,p. 



36pr' + 3(1? • V's) 



+i,i-i-y'sr-l(v'l) 



(14a) 



the Navier-Stokes equation, 

-l^Ps(y''s - V\s+l)mod2) + Fq ■ 

the state equation, 

Ps = isp: , (21) 

where Ps is the fluid pressure, and the Maxwell equations. 
The second ones are governed by the state equation H21(l , 
Maxwell equations, the continuity equation 



V • K' = 



(22) 



/o^(^)°^(x,t) = 6p,c2 



i-4^(46pr' + ^'') 



(14b) 



where the weights Wi are wq — giy, Wi,2.3,4 = ji^i 
W5,6 = 3Bc^- I^ addition, ^^ is a constant that is fixed by 
the initial fluid temperature and density by means of the 
ideal gas law. 
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with polytropic index 7, and k is the Boltzmann constant. 
The equilibrium velocity V's is defined by 



v. 



Vs 



XsTsStF'-'^ 
Ps 



For the electromagnetic field {K = 2), we have 
rf(2)oq.^ t) = —E' -^ +-B -if 



(16) 



(17a) 



and the Navier-Stokes equation for an incompressible and 
viscous fiuid. 



Ps[^ + iV'-V)V'^ 



\ dt 



VK 



+ ^PsiE + V'sXB) (23) 



- lypsiV's - V' (s+l},nod2) 
+ Fo+T]sPsy^V's . 

where the kinematic viscosity is ?7s = f (ts — l/2)c'^St. 

III. SIMULATION OF A 2D HARTMANN FLOW 

In the MHD limit, the two- fluid theory becomes the 
MHD (one fluid) theory, which is represented by the fol- 
lowing equations: the continuity of mass, 



v.(py) + | = o , 



(24) 



p(2)cq 



(17b) 



(18a) 



where the equilibrium electric held is 

E' = E-{noc^X2T2St)J' , 

and X2 — ^, as before. 

The proof that this lattice Boltzmann model, via a 
Chapman-Enskog expansion, recovers the equations of 
the two-fluids theory for a plasma composed by electrons 
and ions is shown in Appendix ^ The model let us 
to consider either compressible and non- viscous fluids or 



the Navier-Stokes equation. 



d 



p{^ + V -V ]V = -VP + J X B + 'qV^V + Fo., (25) 



the magnetic field equation, 
dB 



dt 



\/ X {V X B)+ rim^^B 



and the state equation, 



(26) 



(27) 



where p is the total mass density, V is the total velocity 
field and 77™= — ^ — is the magnetic viscosity. 

For the Hartmann flow [la . Ilq |. we consider a fluid 
in isotermal equilibrium (7 = 1) at low temperature (a 
small ^s value), incompressible and viscous. The fluid 
moves in the x direction between two walls at rest at 
y=—L and y=—L. There is a constant magnetic field 
in the y direction, with intensity Bq, and a constant ex- 
ternal force F—pg in the x direction to drag the fluid 
[l5|. So, the velocity and magnetic fields take the forms 
y=(K(2/),0,0) and B = (B,(y), Bq, 0), respectively. By 
replacing these expressions in equations l|^ and (PS|l . 
one finds the following solutions for the velocity and mag- 
netic fields [l5l: 



V^{y) 
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cosh(iJ) 
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smh{H) L 



,(28a) 



(28b) 



- BqL 



is the Hartmann number and ~L < y < 



where H= 
L. 

For the simulation, we use a single row of 80 cells in 
the y direction, with periodic boundary conditions in 
both X and z directions. The initial conditions for the 
density functions are obtained from the equilibrium ex- 
pressions (I14II and (|17|) with the values 14=0, ps=msns, 
i5=0, B={0,Bo,0) and ^o^"^ = (p^5,0, 0). In addition, the 
constant values are 7=1, Cs=3 x 10~^, /xo=1.0, c=l, 
1^=100, T^=1.0, T2=0.5, mo=1.0 X IQ-^^, mi = 1820mo, 
and no=ni=1.0 x 10^^ particles per unit volume. For 
the y direction, we assume as boundary conditions at 
the walls that the equilibrium density functions for the 
time evolution (Eq. H14|) and H17|) ') are always the same 
from the initial conditions (including T4=0, i.e. non- 
conducting walls). The system evolves until a steady 
state is reached. We ran simulations for Hartmann num- 
bers H=5, 13 and 26, and the magnetic field Bo was 
chosen to obtain these Hartmann numbers. 

Figure [SI shows the velocity profiles and figure El shows 
the magnetic field profiles for the three cases. The solid 
lines are the analytic solutions fEa. H28|) '). The simulation 
results are in excellent agreement with the analytical so- 
lutions. This result say us that (at least for the MHD 
limit) our LB models works properly. 



IV. APPLICATION TO MAGNETIC 
RECONNECTION 

A. Dynamics of the magnetic reconnection process 

In order to simulate the magnetic reconnection in the 
magnetotail, we chose the initial equilibrium condition 
proposed by Harris [I3, [Q for the current sheet, plus 
a magnetic dipole field, ortogonal to the sheet. For this 




FIG. 5: Velocity profile Vx vs. y/L for difi'erent Hartmann 
numbers; H=6.0 (circles), H=13.0 (squares) and H=26.0 (di- 
amonds). The solid lines are the analytical results. 
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FIG. 6: Magnetic field intensity Bx vs. y/L for different 
Hartmann numbers: H=6.0 (circles), H=13.0 (squares) and 
H=26.0 (diamonds). The solid lines are the analytical results. 



simulation we assume that the fluids are non- viscous and 
compressible. 

The current sheet lies on the x-y plane, and its 
magnetic field is described by the vector potential 

A={0,Ay,0), with 



Ay{x,z) = LBol'acosh[v{x){z/L)]/v{x) 



(29) 



where the effective thickness of the current sheet is given 
by L/v{x), and the asymptotic strength, Bq, is the value 
of Bx in the limit z —^ 00, divided by v{x). The function 



v{x) is an arbitrary slowly-varying function. We choose 
for v(x) the quasi-parabolic function proposed by pjll2Clj| . 



v{x) = exp(— ex/L) , 



(30) 



where the parameter e is much smaller than one and de- 
termines the strength of the z-component of the magnetic 
field. We took e=0.1 for the simulation. The initial den- 
sity is the one proposed by Harris, 

ns{x, z) = Ub + ncv'^{x) cosh^ [v{x){z/L)] , (31) 

where Ub is the background density and Ub + ric is the 
maximal density. 

The magnetic dipole is set at position xq with momen- 
tum M and oriented in the z direction. It generates a 
magnetic field given by 



Bx{x,z) = 



3M{x — xa)z 
{{x - xoy + z'^)i 



By{x, z) =0 



(32) 



B,{x, 



M{2z^ ~{x- xof) 
{{x - xo)^ + z^)i 



The lattice constant Sx is chosen as one seventh of 
the ion inertial length, Sx—jc/uji, where wi is the ion 



plasma frequency, Wi 



_ /gfru 



, with ni=10'^ particles per 



cubic meter for the magnetotail |21l | and mi the proton 
mass. That gives (5a;~103km. Since the current sheet in 
the magnetotail can be assumed around 3000km width 
[2l|,|22|, we chose L=2c/uji. For the position of the mag- 
netic dipole, we took xq=22.1c/loi and for the dipole mo- 
mentum, A/=3 X 10^^. The grid is an array of 100x100 
cells on the x-z plane with periodic boundary conditions 
in the y direction and free boundary conditions for the 
fields in the other directions (each boundary cell copies 
the density functions of its first neighbohr in ortogonal 
direction to the boundary at each time step). Thus, the 
simulation region is a square of 14.26c/wi length (around 
10300km). For this simulation we took too=?71i/100 (i.e. 
an electron mass 20 times larger than the real one) in or- 
der to obtain numerical stability, but it has been shown 
|2^ that this point does not qualitatively change the 
physical results. The temperature ratio is chosen to be 
'?o/Ti=0.2, acording to observational results |2J]. For 
this simulation, we took nc=5nb and n(,=0.17ni. Fig- 
ures 13 IHl El and ^1 show the evolution of the magnetic 
field lines in the magnetic reconnection process. This 
appears in a natural way, without the a priori introduc- 
tion of any resistive region. The factor fii is the ionic 
cyclotron frequency, fii = qiBo/mi. This result tell us 
that the model can actually simulate the magnetic recon- 
nection. This simulation took Ih in a Pentium IV PC of 
2.8GHz, i.e. it is really fast. 
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FIG. 7: Magnetic field lines in the magnetic reconnection 
process at t=0 (initial conditions) 
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FIG. 8: Evolution for the Magnetic field lines in the magnetic 
reconnection process, at f = S/fii 
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FIG. 9: Evolution for the Magnetic field lines in the magnetic 
reconnection process, at f = 15/ili 
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FIG. 10: Evolution for the Magnetic field lines in the magnetic 
reconnection process, at f = 20/r2i 
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FIG. 11: Magnetic field lines in the magnetic reconnection 
process at t=0 (initial conditions) 



FIG. 12: Evolution for the Magnetic field lines in the magnetic 
reconnection process, at f = 1.92ms 



B. Reconnection rates 

To compute real reconnection rates we performed a 
similar simulation to the one before, but with the ac- 
tual ratio between electronic and ionic masses (mi = 
1820too)- This choice bring us to take a shorter time 
steps {St~3.76 x lO^^s) and smaller cells ((5x=15.95km) 
in order to reproduce with accuracy the electron moves. 
The LB array is 200x100 cells (larger in direction x), 
for a total simulation region of 3190km in x and 1595km 
in z. Since the region is smaller than before, v{x)—l is a 
good approximation on the entire region. The simulation 
constants are L=1595km ^ and Bo=10.0nT [pL The 
densities in Ea. (|^ are nb=0 and nc=10^m~'^ |23, the 
electronic temperature is chosen as To=5.8MK and the 
ionic one as Ti—23.2AdK |2^. All these are observational 
data. The electronic mass is taken mo=9.11 x lO^'^^kg 
and the ionic mass is toi=1.67 x lO^^^kg. All other con- 
stants of our LB model take their standard values in IS 
units. 

The initial configuration of the magnetic field is shown 
in figure [TTl and the same field after t=1.92ms is shown 
in figure [T^ The reconnection rate we obtain from this 
simulation is i?=0.109, which is in good agreement with 
the experimental observations around i?~0.1 |23. This 
simulation took just 5 minutes in a Pentium IV PC of 
2.8GHz. 



V. CONCLUSION 

In this paper we introduce a 3D lattice Boltzmann 
model for simulating plasmas, which is able to simulate 
magnetic reconnection without any previous assumption 
of a resistive region or an anomalous resistivity. The 
model simulates the plasma as two fluids (one electronic 
and one ionic) with an interaction term, and reproduces 



in the continuous limit the equations of the two-fluids 
theory and, therefore, the MHD-Hall equations. This 
model can simulate either conducting and viscous flu- 
ids in the incompressible limit or non-viscous compress- 
ible fluids, and sucessfuUy reproduces both the Hart- 
mann flow and the magnetic reconnection in the magne- 
totail. The reconnection rate we obtain with this model 
is i?=0.109, which is in excellent agreement with obser- 
vations. 

Since this method includes both electric and magnetic 
fields, plus the density and velocity fields for each fluid, it 
gives much more information on the details of the plasma 
physics. Moreover, this opens the door to much more so- 
phisticated boundary conditions, like conductive walls or 
electromagnetic waves in plasmas. This is an advantage 
upon other magnetohydrodynamic LB models. Further- 
more, it is 3D, so many interest phenomena can be inves- 
tigated here. The model does not require large compu- 
tational resorces. It just takes between 5 minutes and Ih 
in a Pentium IV PC of 2.8GHz and uses around 100MB 
of RAM. 

The model introduces the forces at first order in time, 
but this is not a problem for weak electromagnetic fields 
and low resistive plasmas. If this is not the case, it is pos- 
sible to modify the charge/mass ratio, but this changes 
the MHD-Hall equations and slows the evolution of the 
electromagnetic fields. Another way to increase the nu- 
merical stability consists of modifying the model to re- 
produce the two fluids in a different way: by defining 
density functions for the sum, f^ + f^ , and the dif- 
ference, *L/P(o) + ^yP(i) of the two fluids. It is also 
possible to develop a LB model with 13 velocity vectors 
for the fluids, as proposed by [23. These are promisory 
paths of future work. 

Hereby we have introduced a 3D lattice Botzmann 
model that reproduces the two-fluid theory and includes 
in a natural way many aspects of interest in plasma 



physics, like electric fields and magnetic reconnection. It 
has been shown in this work that this model can actually 
be used to investigate real astrophysical problems. We 
hope that this LB model will contribute to the study of 
plasma physics in many interesting phenomena. 
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APPENDIX A: CHAPMAN-ENSKOG 
EXPANSION 

The Boltzmann equations for each fluid, Eq. ©, © 
and H10|l . determine the system evolution. This evolution 
rule gives in the continuum limit the macroscopic differ- 
ential equation that the system satisfies. This is known as 
the Chapman-Enskog expansion. To develop it, we start 
by taking the Taylor expansion of these equations until 
second order in spatial and temporal {6x, 6t) variables. 
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It is assumed that only the 0th order terms in e of the 
distribution functions contribute to the macroscopic vari- 
ables. So, for n > we have 
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The external forces F*-*-* and the current density J' are 

(Al) of order e [lj|) so we can write F^'^'^—eF^''' and J'=eJ'i. 

Because fp'^'^'^ and flj are now functions of F'^^^ and 

J', we need to develop a Chapman-Enskog expansion of 
the equilibrium function, too: 
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Thus, by replacing these results into Eas. ljAip . (|A2|I 
and HA3|I . we obtain at zeroth order of e 
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where a, /3=a;, y, z denotes the components in x, y and 
z directions. 

Next, we expand the distribution functions and the 
spatial and time derivatives in a power series on a small 
parameter, e, 

^(2)^^(2)(0)^^^(2)(l)_^^2^(2)(2)^^^^ (^4) 
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For the first order terms in e of the distribution func- 
tions we obtain 
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From these equations we can obtain 
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and for the second order terms in e we have 



(Al2c) 






K^St 



d 



i_/j-P(s)(2) _ rf(s)(2)cq^ 



(Al3a) 



9A,r, 



+ 



^{^^.VsM-f['^) (A15b) 



2c2 



3A,r,(F, -i^i^ 



(«)^ 



/, 



'P(s)(2)oq/-; 



(a-,t) 



CJ.f^t 



_^\2^2^,-^ p('^)^2 

4c2' 



A^,^(^f.F^)^ 



'^A^,^(^(^V 



(Al5c) 



and 
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The terms of order one and two for the cquihbrium 
functions of the fluids are obtained by replacing Eq. I|16|) 
into Eq.(^. That gives 
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/o^(^)(")'='^(f,i) = 



6p,c2(l--ij(46pr^ + (V^.)2)^ ^^^^"^^ 



/o^(^)(^)'='^(a-,i) = ^6Stc^^{V,-Ft^)], (Al5e) 



/oP(^)(^)°^(x,t) = 






Ps V 4c2 



?(s) p(s)~ 



(A15f) 



The same process can be used to determine the terms 
of order one and two for the equihbrium functions of the 
electromagnetic fields. Replacing Eq. (|18|l into Eq. p7|) 
and grouping, we have 



/f)(°)<=^(f,t) = ^^.er,. + is.6f^. , (A16a) 
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Now, we are ready to determine the equation that the 
model satisfies in the continuum limit. First, let us con- 
sider non-viscous compressible fluids, that is Ts=^. By 



summing up Eq. (jA12ap over i and p, and by taking into 
account Eqs. (|A12c|l . Q, (|A15p and l|X8)l . we get 
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By summing up Eq. (|A13a|l in the same way, we obtain 

'Xs + l^s j.,^{s)\ , dps 
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Now, we can add these two equations to obtain 
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■-StF, 



. (A19) 



Next, foUowing Buick and Created [Tj|, we do X^ — ^, 
Ats= ^2°T~^ and, by taking into account Eq. (|16|l . we arrive 
to the continuity equation 
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By multiplying Eq. HA12al) by vf and summing up 
over i and p, we get 

^(P.F,„V;,) + ^— + ^^— - F,„ . (A21) 



In a similar way, by multiplying Eq. ljA13a|) by vf and 
summing up over i and p, we obtain 
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Now, we can add these two equations, and by replacing 
Eq. ((TCll . we get (up to second order in e) 



dip.V^J , d 



dP, 



This is the Navier-Stokes equation for non-viscous com- 
pressible fluids, with state equation Ps=^sP2- In our 
model, the force Fa is taken at first order in e. With this 
approximation, Ea. lfTT|l gives F^'^J {Vs)=fI^J {V' s), and 
the Navier-Stokes equation is 



dP / n ^ ^ 



VPsiV's - V\s+l)mod2) + Foa 



By replacing Ea. (|A20p into Ea. (|A24p . we arrive to the 
usual form of the Navier-Stokes equation for a non- 
viscous compressible fluid |j] 



\ dt 
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rus 
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Second, let us consider both fluids with viscosity (t, > 
1/2) in the incompressible limit. By following the same 
procedure, we arrive to the following momentum equa- 
tion (up to second order in e): 

djPsVL) , d 

-^+(^PsiE + V'sxB) ^^ ^ 
dxa \ms (A26) 

- t^Ps(V''s - V\s+l)mod2) 

+ VsPs^^V^c + Foc , 



where the kinematic viscosity is r]s—^{Ts — l/2)c^5t. By 
following the same procedure described above Q, we ar- 
rive 

-~^Ps + ^ME + V'sXB) (A27) 

- I^PsiV's -V\s+l)mod2) 
+ Fq+7],PsV^V's . 

For the electromagnetic fleld, we take T2 = 1/2, A2 = 1 and 
K2 — 0. By summing up Eqs. (|A12bl) and (|A13b|) on 
i, j and p, we do not get any information about the 
fields. Thus, let us multiply these equations by ef before 
summing up. So, we obtain 

dE 
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c^V X B = -p.oc^J'i , 
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If we add these two equations, and because of Eq. l(TF)>. 
we get the first Maxwell equation. 
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Similarly, multiplying Eqs. ljA12b|l and (jA13b|) by fef^ 
and summing up on i, j and p, we obtain 
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If we add these two equations, we obtain the second 
Maxwell equation. 
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The other two MaxweU equations can be obtained from 
the Eas. HA30() and (|A33() as foUows ^. If one apphcs the 
divergence to these equations we obtain 
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Now, we replace the Eq.|(T21) in the Ea. l|A34(l to get 
d{V ■ E') _ 



dt 
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and because of the two fluids satisfy the continuity equa- 
tions IIA20II. we obtain 
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By taking into account the Eq. Q, we finally get 
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Thus, if the initial conditions for the electromagnetic 
fields satisfy the Maxwell equations 



V-B = 



V ■ E' 



2 Pc 

/ioC Pc = — 

eo 



(A39) 



(A40) 



this equations will be recovered for all times. 

Summarizing, the state equation and Eqs. IJA20|I . 
(|A24|I determine the behavior of a non- viscous compress- 
ible plasma. If we use Ea. ljA26p instead of Ea. (|A24|l . 
the model reproduces the behavior of an incompress- 
ible plasma with viscosity. Eqs. (|A30|I . HA33P (|A39|) 
and (|A40p determine the evolution of the electromagnetic 
fields. These are the equations of the two-fluids theory 
(^], and this completes the proof. 
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